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Abstract 

A unified picture of nucleon antinucleon annihilation into pions 
emerges from a classical description of the pion wave produced in an- 
nihilation and the subsequent quantization of that wave as a coherent 
state. When the constraints of energy-momentum and iso-spin con- 
servation are imposed on the coherent state, the pion number distri- 
bution and charge ratios are found to be in excellent agreement with 
experiment. 
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Nucleoli antinucleon annihilation at rest or at low energy goes mostly into 
pions. The problem for theory is to calculate the pion spectrum, the average 
number of pions, the branching rate into various pion modes, higher moments 
of the pion distribution, correlations and charge ratios, all from some theory 
of the annihilation process. In this note we show that a unified account of all 
these features comes naturally out of a picture of very rapid annihilation into 
a classical pion wave and subsequent quantum description of this wave as a 
coherent state. A very simple parameterization of the annihilation process 
coupled with the constraints of iso-spin and four-momentum conservation 
then leads to an excellent account of the principal features of the annihilation 
spectrum. 

Recent studies of annihilation in the Skyrme model H have suggested 
that annihilation proceeds very rapidly when the baryon and anti-baryon 
collide and that the product of this rapid annihilation is a pion pulse or 
coherent pion wave. This is a classical picture. We can quantize that wave by 
using the method of coherent states ||. The relatively large energy released 
in annihilation at rest (13.9 pion masses) and the moderately large number 
of pions seen on average (~ 5) makes a classical starting point plausible. The 
empirical success of the coherent state description gives the quantized rapid 
pion pulse account credibility. 

Suppose the pion wave from annihilation has a source S(f, t), that is well 
localized in space and time, with a corresponding Fourier transform s(k,u). 
If we assume the pion wave obeys a linear wave equation with this source, 
the pion radiation field, <j)(r, t) for t > is given by 

<j>(r, t) = -J ^-s(k, u k )e^-^ (1) 

where uj k = \/k 2 + fi 2 and fi is the pion mass. Introducing creation and 
annihilation operators for each mode k, at and ag, one can define the coherent 
quantum state corresponding to <f)(r, t) by 

|/) = e -^ 2 e/ /( ^ 3fc 4|0) (2) 

with 

fit) = -fc/S^. (3) 
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The mean number of pions in the state |/) is 

N = f\f(k)\ 2 d 3 k (4) 

and|/(/c)| 2 is easily seen to be the single pion momentum distribution in that 
state. The mean energy released is 

E = J u k \f(k)\ 2 d 3 k (5) 

We will set this energy equal to the energy released in annihilation. 

The coherent state contains all possible numbers of pions and hence treats 
all pion annihilation channels together and in a unified way. These channels 
are distributed in a Poisson distribution with mean number N and variance 
a = We arrange parameters so that N ~ 6. We have chosen 6 rather 

than the experimental 5 to roughly compensate for decay into resonant chan- 
nels. These contribute to making N small, but are not part of our "pions 
only" picture. With this N, we find too large a variance. The experimen- 
tal result is u ~ 1. Furthermore only states with between 2 and 13 pions 
are permitted in annihilation at rest by energy and momentum conservation, 
while the coherent state contains any number of pions. To maintain the si- 
multaneous treatment of all annihilation channels that is characteristic of the 
coherent states, and correct its failings, we need to introduce the constraints 
of energy and momentum |4j]. To do this, introduce the operator F(x), 

F(x) = J d 3 kf{k)ate ik - x (6) 

that creates a state at the four-vector position x with Fourier components 
f(k) (the same f(k) as above). In (6) the time component of k is the on 
shell energy u^. The state of fixed total four- momentum K is then given by 

\f,K)= f ^e- tK - X e F(x) \U) (7) 

as is easily checked by expanding out the exponential e F ( x \ One finds 
(/, K'\f, K) = 5\K - K')I(K), with 

I{K) = J T^^e^) (8) 
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and with 

p(x) = I d"p\f(p)\ 2 e-^ x (9) 



The integral for I(K) is singular because of the large x behavior of the factor 
e p( x )_ This singular behavior comes from the "one" in the expansion of this 
exponent which comes in turn from the "one" in the expansion of e F<yX \ This 
is the no pion state and is forbidden by four-momentum conservation, as is 
the one pion state. Hence we can write 



(2vr) 4 ^ 2 m\ 

= E^r (10) 

The sum over m terminates with m = 13 by four-momentum conservation. 
The I m (K) are easily calculated by first doing the x integration, extracting 
the four-momentum conserving delta function and calculating the remaining 
integrals with the constraint imposed. In terms of the I m , the probability of 
finding m pions is P m = jf^^\ and the mean number of pions is given by 

n = y mIm{K) (in 

To calculate P m and N, we need to model S(f,t). Inspired by the 
Skyrmion calculations, [|TJ we take a very simple spherically symmetric 
form. We assume S(f, t) = for t < and for t > 0, 

S(r, t) = Cte-^e- ar /r (12) 

where C is a scale constant. This leads to 

C'k 2 

l/(fc)|2 = (pT^FR+TM (13) 

where C is a another scale constant and where we have multiplied / by k 
to model the p-wave nature of pion emission. We fix C by requiring that 
the average energy be the energy emitted in annihilation at rest, (5), which 
is 13.87 in units of the pion mass, (fi = 1). For the range parameters we 
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make the very simple assumption that = 7, and after very little parameter 
searching find that a = 7 = 2 gives N = 6 for the Poisson distribution. Thus 
this very simple form with the reasonable size of half the pion Compton 
wave length begins to look like the data. Using the same f(k) in the four- 
momentum restricted calculation we find N = 6.4 and a = .88, both in quite 
close agreement with experiment roughly corrected for final state resonances 
1, 1- 

In Figure 1 we show the probability of finding m pions in annihilation at 
rest calculated with the unconstrained form (Poisson distribution), and with 
the constraint of four-momentum conservation. We see that the constraint 
sharply narrows the distribution. In fact on the scale of the figure, the 
constrained distribution is indistinguishable from a Gaussian distribution 
of the same average and variance. The constrained distribution plotted in 
Figure 1 is nearly identical to the distribution seen in experiment 0,0. 

So far we have neglected the iso-spin of the pion. Since for averages, the 
constraint of four- momentum conservation does not seem very important, let 
us examine the effects of iso-spin on pion averages in the ordinary coherent 
state. We must now make the pion creation operator an iso-vector. But we 
cannot make f(k) an iso-vector and dot it into the creation operator because 
iso-spin conservation is a global constraint. It must hold for every k. To 
implement iso-spin conservation, we introduce a fixed unit vector in iso-spin 
space, T, dotted into the creation operator, and write the coherent state with 
fixed total iso-spin, J, and z-component I z as 

\f,I,I.) = v I ^Y^(f)e* J '^-l^O) (14) 

where v is a normalization factor, and N is the average number of pions 
summed over charge types. It is still given by (4). The expectation value of 
the number operator for pions of charge type \x in this state gives the average 
number of pions of this type. For 1 = one easily finds the expected answer 
of N^ = N/3, for any fi. For I = 1, I z = on finds, N + = N_ = N/5 and 
Nq = 3N/5. These results are obtained in the large N limit where they are 
independent of the form of f(k). Note the excess of 7To in the 1 = 1 case. 
If we use a recent estimate of the relative population of I = and / = 1 in 
proton antiproton annihilation at rest ||, we find Nq/N + = 1.53 ± .15 (the 
uncertainty comes from the population estimate), to be compared with the 
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experimental value of 1.27 ± .14 0,0. Thus the coherent state constrained 
by iso-spin conservation agrees remarkably well with the data for the pion 
charge ratios and naturally accounts for the excess of 7r 's. 

We have seen that a description of proton antiproton annihilation at rest 
into a classical pion wave and the subsequent quantization of that wave as 
a coherent state gives a unified picture of all the direct pion channels and 
correctly accounts for the number distributions of those channels and for the 
charge ratios, all with very few parameters. Imposing four-momentum and 
iso-spin conservation on the coherent state is important for this agreement. 
These constraints also introduce correlations that are not in the simple co- 
herent state. Extensions of these methods to study these correlations and 
to other conservation laws are under study. Further afield one can imagine 
applying these ideas to any strongly interacting process in which the origi- 
nal event might be described classically and quantum mechanics restored by 
coherent states. Jets, Centauros, and very high energy heavy ion collisions 
may be three such cases. For Centauros and heavy ions there has recently 
been some work along these lines ||, [[UJ, [|TTJ . These approaches, like ours, 



take classical, nonperturbative QCD as their starting point. Much of this 
work uses polynomial coherent states with fixed numbers of pions in order to 
introduce quantum mechanics and to implement iso-spin conservation. The 
very large number of pions (> 100) seen in these processes suggests that 
field coherent states with our prescription for iso-spin are a more natural 
approach. One might imagine many other processes in which QCD can be 
treated classically and then experimental manifestations described by coher- 
ent states. In many cases, including annihilation, it may be necessary to 
further develop the treatment either by introducing squeezed states or by 
using a density matrix to average over coherent states when processes have 
important statistical components as in the case of heavy ions. 

RDA, FC, and JPD would like to thank the theory group of the Division 
of Nuclear and Particle Physics at the Paul Scherrer Institute for providing 
a very pleasant environment in which much of this work was done. RDA and 
BS are supported in part by the United States National Science Foundation. 
The Division de Physique Theorique is a Research Unit of the Universities 
Paris 11 and 6 associated to CNRS. 



6 



References 

[1] H.M. Sommermann, R. Seki, S. Larson and S.E. Koonin, Phys.Rev. 
D45, 4303 (1992). 

[2] B.Shao, N.R.Walet, and R.D.Amado, Phys.Lett. B 303, 1 (1993). 

[3] R.J. Glauber, Phys. Rev. 131, 2767(1963). 

[4] D. Horn and R. Silver, Annals of Physics 66, 509 (1971). 

[5] J. Sedlak and V. Simak, Sov. J. Part. Nucl. 19, 191 (1988). 

[6] C.B. Dover, T. Gutsche, M. Maruyama and A. Faessler, Prog. Part. 
Nucl. Phys. 29, 87 (1992). 

[7] C. Ghesquiere in: Proc.Int.Symp. on Antinucleon-Nucleon Interaction, 
Liblice-Prague, 1974, CERN 74-18, p. 436 

[8] M.P Locher and B.S. Zou, Z.Phys A 341, 191 (1992). 

[9] C. Greiner, C. Gong and B. Muller, Phys.Lett. B 316, 226 (1993). 

[10] J-P. Blaizot and A. Krzywicki, Phys. Rev. D46, 246 (1992). 

[11] K. Rajagopal and F. Wilczek, Nucl. Phys. B404, 577 (1993). 



Figure Caption 



Fig.l The probability, P m , of having m pions in nucleon antinucleon 
annihilation at rest as a function of m. The open circles refer to the 
unconstrained coherent state and are a Poisson distribution. The solid 
squares include the constraint of four-momentum conservation. This set of 
points agrees well with experiment, 0, 0, and is indistinguishable (on the 
scale of the figure) from a Gaussian distribution of the same mean and 
variance. 
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